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1.  Introduction 


For  many  years,  the  Advanced  Munitions  Concepts  Branch  of  the  U.S.  Army  Research  Labora¬ 
tory’s  (ARL’s)  Weapons  and  Materials  Research  Directorate  has  been  producing  trajectory  heading 
histories  for  flight  experiments  of  developmental  projectiles  using  data  telemetered  from  strap- 
down  optical  and  magnetic  sensors  (Hepner  &  Harkins,  2001).  U.S.  Patent  638155,  Projectile 
Orientation  In  Navigation  TERms  (POINTER),  was  awarded  in  June  2002  for  the  design  of  a  con¬ 
stellation  of  optical  and  magnetic  devices  and  a  processing  methodology  that  provides  a  solution 
for  the  projectile  pitch  and  yaw  angular  history.  The  POINTER  solution  is  obtained  from  the 
simultaneous  solution  of  a  system  of  equations  describing  respective  sensor  output  as  functions  of 
projectile  orientation.  Because  one  of  these  equations  is  quadratic,  there  are  two  possible  solutions. 
This  ambiguity  usually  is  easily  resolvable  by  the  analyst  when  s/he  is  post- flight  processing  the 
data.  However,  automated  and  on-board  processing  sometimes  cannot  identify  the  correct  solution. 
An  alternate  set  of  simultaneous  equations  that  yields  a  unique  POINTER  solution  has  been  identi¬ 
fied  and  is  presented  herein. 

The  “legacy”  POINTER  solution  is  reviewed  in  sections  1  and  2  before  we  proceed  to  the  develop¬ 
ment  of  the  new  POINTER  solution  in  sections  3  and  4. 


2.  Solving  for  Projectile  Heading:  Generalized  Solution 


The  center  of  motion  and  the  principal  axis  of  rotation,  i.e.,  the  spin  axis,  of  a  three-dimensional 
freely  flying  solid  body  are  described  within  a  Cartesian  coordinate  system  with  the  use  of  three 
variables  denoting  location  and  two  variables  defining  angular  orientation.  The  derivative  of  the 
location  with  respect  to  time  is  commonly  referred  to  as  the  velocity  vector,  v  .  Eor  a  symmetric, 
spinning  body,  the  navigation  pointing  vector,  p  ,  normally  is  coincidental  with  the  principal  axis 
of  rotation.  The  rotation  rate  of  the  body  about  this  axis  is  commonly  called  the  “spin  rate”.  The 
Eulerian  heading  variables  psi  (i//)  and  theta  (6)  are  used  to  denote  the  two  angular  components  of 
azimuth  and  elevation  required  to  orient  the  body  principal  axis  of  rotation  within  a  right-handed 
Cartesian  navigation  system.  The  principal  axis  of  rotation  of  a  spinning  body  is  often  not  co- 
linear  with  the  velocity  vector.  In  such  cases,  an  orientation  time  history  estimated  from  deriva¬ 
tives  of  location  variables  does  not  provide  an  accurate  measure  of  the  navigation  pointing  angles 
of  the  flight  body.  With  the  combination  of  body-fixed  sensors’  measurements  of  the  included 
angles  between  the  axis  of  rotation  and  two  distinct  earth-fixed  fields  and  knowledge  of  these 
fields’  orientations  in  the  navigation  coordinate  system,  ^and  0  can  be  determined.  Thus,  a 
pointing  angle  time  history  can  be  generated. 
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The  parameters  just  described  are  all  seen  in  figure  1  where  a  projectile  with  its  pointing  and 
velocity  vectors  is  shown  in  an  exemplary  right-handed  coordinate  system  along  with  two 
representative  field  vectors  of  known  orientation,  ,  and  .  The  included  angles  between  the 

spin  axis  and  the  two  field  vectors  (ci  and  02)  and  the  azimuth  (i//)  and  elevation  (0)  heading 
angles  are  also  shown. 


Let  the  unit  vectors  P  ,  F^,  and  along  P ,  F^  ,  and  F^  be  defined  within  the  navigation  {X,  Y, 
Z)  system  as 

P={P„Py,P,) 

F,={F,^,F,^,F,^)  (1) 

F,  =  (F,  ,  F,  ,  F,  ) 

The  components  of  P  are  obtained  from  the  simultaneous  solution  of  the  system 

P  •  F^=  cos(crj ) 

P  •  F^  =  cos(cr2 )  (2) 

|f|  =  i 

_  _  -> 

where  the  field  vectors  (Fj ,  F^)  and  the  included  angles  (cri,  02)  between  P  and  the  respective 

field  vectors  are  known,  estimated,  or  measured.  This  methodology  is  termed  POINTER  (U.S. 

Patent  6398155,  June  4,  2002,  Harkins  &  Hepner)  and  has  been  successfully  employed  on 
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telemetered  data  from  numerous  inventory  and  developmental  projectiles  instrumented  with  solar 
and  magnetic  sensors  in  flight  experiments  conducted  to  characterize  projectile  dynamics  (Davis  et 
ah,  2004;  Harkins,  2003). 


3.  Solving  for  Projectile  Heading:  The  Legacy  POINTER  Solution 


Angular  navigation  parameters  are  most  often  described  in  an  earth-fixed  Cartesian  coordinate 
system.  Given  a  convenient  navigation  coordinate  system  (e.g.,  north,  east,  down)  and  an  arbitrary 
location  L  =  (Ln,  Le,  Ld),  let  the  unit  vector  from  L  to  the  sun  be  5  =  {S„  and  the  unit 

vector  from  L  along  the  magnetic  field  be  M  =  ) .  Unfortunately,  the  angular  vari¬ 

ables  in  the  system  of  equations  2  describing  heading  usually  are  non-separable  in  this  coordinate 
system.  In  order  to  obtain  the  projectile  heading  angles,  the  vector  along  the  projectile  axis  must  be 
derived  in  a  computationally  convenient  coordinate  system  and  then  transformed  into  the  navigation 
coordinate  system. 

Consider  a  Cartesian  coordinate  system  (A,B,C)  with  its  origin  at  L,  its  +B  axis  along  the  solar 
vector,  its  A  axis  so  that  the  local  magnetic  vector  is  in  the  half-plane  containing  the  B  axis  and  the 
+A  axis,  and  its  +C  axis  pointing  in  the  direction  that  a  right-hand  threaded  screw  advances  when 
its  head  is  rotated  from  +A  to  +B  (figure  2).  This  system  will  be  referred  to  as  the  “solar”  system. 
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With  the  angle  between  S  and  M  designated  as  p,  the  mappings  of  the  axes  of  the  solar  system 
to  the  navigation  system  are  given  by 

(0,1,0),  5,, 5,)  (3) 

(0,0,1),  ^(Mx5)/sinp  (4) 

(1, 0, 0),  ^  5  X  [(m  X  5)/sin  p]  (5) 

where p  =  eos  '(;S  •  m).  Evaluating  the  eross-produets  in  equations  4  and  5  yields  expressions  for 
the  eomponents  of  the  respeetive  mappings 

M  X  5  =  {M,Sd  -MdS^MdSn  -M„Sd ,  M„S, -M,S„ )  (6) 

.•.(0,0,4  -MdSe,MdSn  -M^Sd,MnSe  -M.^J/sinp  (7) 

5  X  (m  X  5)  =  5  X  {M^S,  -M,S^,  M,S„  -M„S„  ) 


(8) 


(1,0, 0)< 


sXm„s, -mx)-sAm,s„ 

S,  {M^S,  -M,Sj-S„  {M„S,  -  M^S„ ), 

{M,S„  -M„S,)-S, {M^S,  -  M,S, ) 


+sA-S„S,M,  -S^SdMd,  ^ 


MAS^^  +s/]-S^S,M^-S,SdM^,  \  /smp 


Md\Sn^  +S,^]-S^SdM^  -S,SdM, 


(9) 


Any  veetor  defined  in  the  “solar”  system,  F,  =  (v^,v,,Vp),  ean  be  transformed  to  its  representation 
in  the  navigation  system,  Vj^  =  (^« ,  ) ,  by 


^h,i 

h,2 

t  \ 

II 

^2,1 

^2,2 

^2,3 

f^3,l 

h,2 

^3,3  y 

K 


(10) 


where 
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h,l  =  (*i„(Se^+s/)-S„S,M, -S„SjMj  j/sinp 
'1,2  =(l^e{s„^+s/)-S„S,M„ -S,SjMj)/smp 
'U  =  (mj(s/+s/)-S„SjM„ -S,SjM,)/smp 

^2,1  “ 

t2,2=Se  (11) 

h,2  =Sd 

h,l  ={MeSd  -MdSe)l^'^np 
^3,2  ~  (  ^ d^n  ~ ^ n^d  )/^l^P 
^3,3  ~  ~ )/^l^  P 

Given  solar  aspect  angle  (gs)  and  magnetic  aspect  angle  (gm)  histories  derived  from  solarsonde 
and  Magsonde  reductions  (Hepner  &  Harkins,  2001),  these  data  should  be  interpolated  onto  a 
common  time  base.  Then,  at  every  time  step,  the  components  of  a  unit  vector  {P)  along  the 
projectile  spin  axis  can  be  readily  determined  in  the  solar  system  (see  figure  3). 


Figure  3.  Components  of  the  projeetile  pointing  veetor  in  the  “solar”  system 


co^CTs  =  {0X^)*{Pa,Pb^Pc)  =  Pb 
coscTf^  =  (sin p, cos p x) •{Pa^Pb^ Pc)  =  Pa  sin P  +  Pb  cos p 
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(12) 

(13) 


by  substitution, 


-  cos  cr^  cos/? 
sin/? 


(14) 


n  _cosa^  - 


Pc=+^-P/-P/ 


=  + 


1  -  cos  CTj  - 


cos  cr^  -  cos  cr^  cos  /? 

sin/?  j 


or 


^  COS^  <Jm  cos  <J^  cos  (7g  cos  /?  +  COS^  (7g  cos^  /?  ^ 


Pc=+Jsin  o-^- 


sin  /? 


(15) 


(16) 


The  sign  ambiguity  of  Pc  is  usually  easily  resolved  with  knowledge  of  the  initial  navigation  orien¬ 
tation  (e.g.,  launeher  orientation). 


Finally,  with  equation  10,  the  eomponents  in  the  navigation  system  of  a  unit  veetor  along  the  pro- 
jeetile  axis  of  rotation  are  given  by 


{Pn,Pe^Pd)  = 


h,\  h,2  h,2, 

h,\  h,2  h,3 

lv^3,l  h,2  h,3 ) 


{Pa,Pb,Pc) 


and  the  eorresponding  heading  angles  by 


0  =  tan 


-1 


p2  ,  p2 

Wpn  jj 


where  -  nll<0  <  njl 


^  =  tan where  -n<\j/<n 


(17) 


(18) 

(19) 


4.  Coordinate  Transformations  With  the  Use  of  Euler  Angles 


Formulation  of  the  inertial  navigation  problem  for  gun-  and  tube-launehed  projeetiles  requires  the 
use  of  multiple  eoordinate  systems.  Trajeetory  time  histories  are  best  deseribed  in  an  earth-fixed 
eoordinate  system  with  its  origin  at  the  launeher.  Of  neeessity,  strap-down  sensor  measurements 
are  made  in  a  flight-body-fixed  eoordinate  system,  and  target  loeations  are  most  naturally  deseribed 
in  another  earth-fixed  system. 

The  first  eoordinate  system  is  right-handed  Cartesian  (/,  J,  K)  with  its  origin  at  the  launeh  site. 

This  will  be  referred  to  as  the  “earth-fixed”  system  and  the  axes  are  defined  by 

•  The  /  and  Jaxes,  whieh  define  a  plane  tangential  to  the  earth’s  surfaee  at  the  origin. 
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•  The  K axis,  which  is  perpendicular  to  the  earth’s  surface  with  positive  downward,  i.e.,  in 
the  direction  of  gravity. 

•  The  /  axis,  which  is  chosen  so  that  the  centerline  of  the  launcher  is  in  the  /-.Siplane. 


Down-range  travel  is  then  measured  along  the  /axis,  deflection  along  the  Jaxis  (positive  to  the 
right  when  one  is  looking  down  range),  and  altitude  along  the  K  axis  (positive  downward)  (see 
figure  4). 


The  second  system  is  convenient  for  aeroballistic  computations  of  rigid  projectiles’  flights  and  for 
describing  the  locations  and  orientations  of  such  projectiles’  components.  This  system  is  right- 
handed  Cartesian  {i,j,  k)  with  its  origin  at  the  center  of  gravity  (c.g.)  of  the  flight  body.  For  rota¬ 
ting  flight  bodies,  the  projectile-fixed  coordinate  system  usually  has  its  i  axis  along  the  projectile 
axis  of  symmetry,  i.e.,  the  spin  axis  (with  positive  in  the  direction  of  travel  at  launch).  The  j  and  k 
axes  are  then  oriented  so  as  to  complete  the  right-handed  orthogonal  system  (figure  4).  Spin  (/?), 
pitch  {q),  and  yaw  (r)  rates  are  measured  about  these  axes.  This  will  be  referred  to  as  the  “body- 
fixed”  system. 

The  third  coordinate  system  {n,e,d)  is  commonly  employed  to  specify  locations  on  or  near  the 
earth’s  surface,  i.e.,  north,  east,  and  down.  This  will  be  referred  to  as  the  “navigation”  system. 

The  earth-fixed  and  body-fixed  coordinate  systems  are  related  through  an  Euler  rotation  sequence 
beginning  with  a  rotation  of  the  earth-fixed  frame  about  the  K  axis  through  the  yaw  angle  y/.  The 
system  is  then  rotated  about  the  new  J’  axis  through  the  pitch  angle  6.  Finally,  the  system  is  ro¬ 
tated  about  the  i  axis  through  the  roll  angle  (j).  The  two  systems  are  related  by  the  direction  cosine 
matrix  (DCM),  lEb,  with  the  subscript  denoting  earth  fixed  to  body  fixed.  This  transfor-mation 
matrix  is 
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T  = 

^Eb 


-Sa 


“*“  ^y/^^  ^ 6^ (ji 


(20) 


where  c.  is  cos(») ,  and  5.  is  sin(») .  Figure  5  shows  both  the  earth-fixed  and  the  body-fixed 

coordinate  systems  and  the  Euler  angle  relations  between  them.  Transformations  between  any 
two  right-handed  Cartesian  systems  can  be  similarly  defined  with  appropriate  values  for  the  Euler 
angles. 


Figure  5.  Earth-fixed  and  body-fixed  systems  and  the  Euler  angle  rotations. 


5.  Solving  for  Projectile  Heading:  A  Unique  POINTER  Solution 


In  section  3,  it  was  shown  that  the  angular  variables  describing  projectile  heading  can  be  made 
independent  and  therefore  separable  in  an  appropriate  coordinate  system.  In  section  4,  it  was 
shown  that  transformation  between  two  similar  Cartesian  systems  is  described  by  the  Euler  rotation 
sequence.  With  the  choice  of  a  suitable  coordinate  system,  these  two  characteristics  are  combined 
to  achieve  a  system  of  equations  whose  simultaneous  solution  is  unique.  This  “computational” 
coordinate  system  is  somewhat  similar  to  the  “solar”  system  of  section  3  and  is  shown  next.  This 
coordinate  system  (Ic,  Jc,  Kc)  has  its  origin  at  the  projectile  c.g.  It  is  oriented  so  that  the  unit  vector 
to  the  sun  {S)  is  located  along  the  -Kc  axis;  a  unit  vector  (M)  along  the  local  magnetic  vector  is  in 
the  half-plane  containing  the  Kc  axis  and  the  +Ic  axis,  and  its  +Jc  axis  is  pointing  in  the  direction 
that  a  right-hand  threaded  screw  advances  when  its  head  is  rotated  from  +  Kc  to  +  Ic  (figure  6). 
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With  the  angle  between  S  and  M  designated  as  p,  the  mappings  of  the  axes  of  the  computational 


system  to  the  navigation  system  are  given  by 

(0,0-l)c^(5„,5„5^)  (21) 

(0,l,0)c^(Mx5)/sinp  (22) 

(l,0,0)c  ^5x((Mx5)/sin/7)  (23) 

where  p  =  cos~^  •  m)  .  Proceeding  as  in  section  3,  the  components  of  the  respective  mappings 

are 

(0,0,l)c^(-5„,-5,-5^)  (24) 

{W)c^{M,Sd  -M,Sd,M,Sd  -M„Sa,MnSd -M,S^)lsmp  (25) 

'm„(s,^+s/)-s„s,m^-s„s^m„  '  / 

(l,0,0)c^  M^(s„^+s/)-S„S,M„-S,S^M„  Lmp  (26) 
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Thus,  any  vector  defined  in  the  “computational”  system,  V(j  =  (v/^  ,v ],  can  be  transformed 
to  its  representation  in  the  navigation  system, 


V  =T  V  = 

'  N  ^CN  '  C 


1,1 


'•1,2 


h,3 


^2,1  ^2,2  ^2,3 


^3,1  ^3,2 


Vr- 


(27) 


3,3  J 


where 


lu  =  ( W  (5/  +  Sf )  -  S,S,M,  -  S,S,M,)/sm  p 

(s:+S,^)-SAM, -SAM,)ls\np 
>21  =(M^Sj-M^S,)lAap 

(28) 

'2j=(3^,S,-3/A)/s>nP 

h,2  ~  ~^e 

h,3  ~  ~^d 

Continuing  as  before,  solar  aspect  angle  (os)  and  magnetic  aspect  angle  (ctm)  histories  derived 
from  solarsonde  and  Magsonde  reductions  should  be  interpolated  onto  a  common  time  base. 
Then,  at  every  time  step,  the  components  of  a  unit  vector  (Pc)  along  the  projectile  spin  axis  can 
be  determined  in  the  computational  system  (see  figure  7). 


As  before,  equations  for  the  solar  and  magnetic  aspect  angles  lead  to  solutions  for  two  of  the 
pointing  vector  components. 


cos  =  (0, 0,-l)  •  ,  Pj^  ,PkJ=  -Pk, 

cos  ctm  =  (sin  p,  0,-  cos  p)  •  ,  Pj  ,  j  =  Pi  sin  p  -  Pjc  cos  p 


(29) 

(30) 


by  substitution. 


cos  O'M  ~  cos  cr^  cos  p 
sin  p 


(31) 


A  solution  for  the  third  component  follows  from  the  equations  for  the  Euler  sequence  necessary 
for  transformation  from  the  computational  system  to  the  body-fixed  system  as  shown  in  figure  8 
The  transformation  matrix  from  the  computational  system  to  the  body-fixed  system  is  given  by 
equation  20,  mutatis  mutandis. 


1 

¥ch  ^Ch 

^¥cb^^ch 

^Ch 

II 

^¥ch^  ^ch  ^¥cb^$ch 

^ ¥cb^ ff>ch  '^^¥cb^^cb 

1 

\^  ¥  Ch^  ^Ch  '^^¥ch^^ch 

^  ¥cb^  ^cb^  $ch  ^¥cb^  ^ch 

^^Cb^^Cb  J 

(32) 
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In  particular,  the  unit  solar  vector,  S  ,  and  the  magnetic  field  vector,  fh^  =M  are  mapped, 

respectively,  into  the  body-fixed  system  by 


(SP 

^¥ch^^cb 

-Sn  ^ 
^Cb 

— 

'  ^¥cb^  ^Cb^  <Pcb  ^¥cb^^Cb 

^  ¥cb  ^  ^Cb  ^  ^Cb  ^  ¥cb  ^  ^Ch 

^^Cb  ^  ^Cb 

0 

ISkJ 

¥ cb^ ^cb^ ^ch  ^¥cb^ <l>cb 

^  ¥cb^  ^cb  ^¥ch^  $cb 

^^Cb^^Cb  ) 

l-lj 

{ 

^¥ch 

^¥cb^^Cb  ^  ^Ch 

Y  sin(/7)  ^ 

m, 


^  ¥  cb^  <l>cb  ^¥cb^4> 


^¥ch^^cb^^ch  ^¥cb^^cb 


\^ ¥ cb^ ^cb^ ^cb  ^¥ch^ (Pcb  ^ ¥cb^ ^cb  ^¥cb^ ^cb 


Cf)  s. 

^Cb  ^ 
C  n  Cl 


0 

-  COS(yC>) 


(33) 


(Kl)  (34) 


By  inspection  of  figures  7  and  8,  it  is  seen  that  Ocb  =  90  -  oj.  For  a  calibrated  tri-axial  magne¬ 
tometer  whose  axes  are  parallel  to  the  body-fixed  system,  equation  34  gives  the  output  of  the 
respective  axes.  The  spin- axis-aligned  magnetometer’s  output  is 


Rearranging  terms  yields 


sin(p)  (|«c|  )+  cos(/7)  [\mc\ ) 

^  mi  -sq^^  cos(/?)(|mc|)^ 

,  sin(p)(|mc|)  ^ 


¥cb  =  cos 


(35) 

(36) 


Because  cos(^)  =  cos(-^),  the  correct  polarity  of  y/cb  must  be  determined  before  we  can  proceed. 
Since  the  optical  detectors  in  the  solar  measurement  system  are  essentially  binary  devices  that  have 
zero  output  when  not  aligned  with  the  sun  and  unit  output  when  aligned,  continuous  values  for  Sj 
and  Sk  are  not  available  f  However,  any  optical  detector  whose  field  of  view  is  coplanar  with  the 
projectile  spin  axis  provides  a  roll  angle  {(/kcb)  index  when  aligned  with  the  sun.  Typically,  there 
are  two  such  detectors  oriented  along  the  +k  and  -k  body  axes.  If  we  consider  the  +A:-oriented 
detector  when  it  is  aligned  with  the  sun,  it  follows  from  equation  33  that  =  180°.  With  equa¬ 
tion  34  for  the  corresponding  values  of  either  mj  or  ink  at  this  time,  the  correct  polarity  of  y/cb  is 
identified.  After  this  polarity  has  been  determined,  it  needs  only  to  be  rechecked  if  and  when 
|cos(^^ci)|  ^  1  • 

With  the  correct  value  of  y/cb  in  hand  and  observing  from  figures  7  and  8  that  tan(^(^^  ]  =  Pj ^  /  Pj^  , 
we  compute 

Pj^  =  tan(^^cfo )  Pp  (37) 


Finally,  with  equations  27  and  28,  the  pointing  vector  in  the  navigation  system  is  uniquely  given  by 


P  =T  P 

^  N  ^CN  ^  C 


(38) 


'if  continuous  tri-axial  measurements  with  respect  to  two  earth-fixed  fields  are  available,  the  possibility  of  com¬ 
puting  Euler  angles  in  multiple  ways  using  different  combinations  of  the  sensor  output  can  be  exploited  to  identify 
and  correct  for  individual  sensor  corruption  and/or  failure  (Jagadish  &  Chang  2007). 
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Elevation  and  azimuth  in  the  navigation  system  are  then  eomputed  with  equations  18  and  19, 
respectively. 


6.  Summary 


A  system  of  equations  relating  strap-down  optical  and  magnetic  sensor  output  and  projectile 
heading  have  been  identified  whose  simultaneous  solution  uniquely  identifies  the  projectile 
azimuth  and  elevation.  This  represents  an  improvement  of  an  earlier  formulation  that  led  to  an 
ambiguous  solution.  With  a  unique  solution  available,  the  possibility  of  on-board  and  automated 
computation  of  projectile  heading  is  opened.  This  should  be  attempted  as  soon  as  practicable  for 
test  and  evaluation  purposes. 

For  the  future,  there  are  a  number  of  potential  alternatives  to  the  solar  field  measurements 
discussed  herein,  which  when  successfully  realized,  will  enable  a  tactical  implementation  of  a 
POINTER  system  using  low-cost  sensors. 
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